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Abstract 

We consider a system made up of exictations of a neutral scalar field, </>, having a Xcf)'^ interaction 
term. Starting from an ensemble where the occupation number / is large, but A/ is small, we 
develop a classical field theory description of the evolution of the system toward equilibrium. A 
Boltzmann equation naturally emerges in this description and we show by explicit calculation that 
the collision term is the same as that coming from elastic scattering. This shows the equivalence 
of a Boltzmann equation description and a classical field theory description of the same system. 
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I. INTRODUCTION 



It is believed that in the very early stages after a heavy ion collision, gluon occupation 
numbers are as large as 1/a, ~ l/a [|l], ^ ^, ^, ^. As the QCD quark-gluon plasma evolves 
toward equilibrium, occupation numbers decrease until ~ 1 at equlibrium. Understanding 
the evolution of dense quark-gluon matter toward equilibrium is both practically important 
and theoretically challenging. So long as <^ l/a the Boltzmann equation, with elastic 
and inelastic collision terms and quantum-statistical factors, should furnish a systematic 
theoretical framework for this problem 0, although this equation cannot be expected to 
be useful at the earliest times after the collision when fg ~ 1/a. On the other hand, when 
occupation numbers are large one can expect classical field theory to apply, and there is 
currently an interesting program studying the early stages after a heavy ion collision using 
classical Yang-Mills equations with an initial condition fixed by the McLerran-Venugopalan 
saturation model [|^, || . Thus one might expect that classical field theory should be a good 
theoretical framework for studying a dense gluon system up to times just before equilibration 
occurs. It would seem then that classical field theory and the Boltzmann equation are 
equivalent frameworks for studying dense non-equilibrium QCD in the period where 1 <^ 
/g ^ 1/a holds for the important phase space region of the system. 

While this equivalence have been studied previously, most notably in connection with 
wave turbulence ||^, it is less familiar in the context of the early stages of heavy ion collision. 
We thus explore this equivalence in detail in this paper. To simplify the discussion we deal 
with a 0*^ field theory rather than with QCD, but we believe our argumentation should 
work equally well for Yang-Mills theories. Our goal is to derive the collision term in the 
Boltzmann equation using only classical field theory, and we do this explicitly at lowest order 
in the coupling for the elastic-scattering part of the collision term. This part of the collision 
term, illustrated in Fig. (|^) below, has a contribution cubic in the occupation number and 
a contribution which is quadratic. Perhaps surprisingly we are able to reproduce both of 
these parts from the classical theory, and the reasons for this are discusssed briefly at the 
end of Sec. ^V\. For inelastic parts of the collision term the classical field description cannot 
be expected to give the complete answer, even at the lowest nontrivial term in the coupling, 
but it should give these portions having the maximum and maximum minus one powers of 
the occupation number. 
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In Sec. m we review the description of a non-equilibrium system, described by a 0^ field 
theory including the usual doubling of the number of fields. When the fields are large we 
show how a classical field theory naturally replaces the quantum field theory description 
of the system, and we introduce combinations of the field variables to make the classical 
description transparent. 

In Sec. jn|we introduce Green's functions in the classical field theory and write equations 
of motions for the Green's functions. If the medium does not vary too rapidly we show how 
a Boltzmann-like equation emerges with a "collision term" given by Green's functions of the 
classical theory, depending of course on the initial ensemble defining the system. 

In Sec. 1^ we identify the "collision" term derived in the classical theory with the usual 
collision term given by elastic scattering. The identification is done by explicit calculation 
at lowest nontrivial order in perturbation theory. 

Our whole discussion is carried out under conditions where the occupation number / is 
large (/ ^ 1) while A/ is small (A/ ^ 1) where A is the usual 0^ coupling constant. We 
believe that higher order corrections in A/ can be done in the classical theory, but higher 
order in A will, in general, be quantum. 



II. DESCRIBING A DENSE NON-EQUILIBRIUM SYSTEM 

The system which concerns us here is made up of excitations of a neutral scalar field 
of mass m and with interactions described by a — interaction term in the Lagrangian 
density. Suppose at time t we wish to determine the expectation of some observable, O, 
made of of field 0. For example, O might be 0^(x, t), (j){xi,t)(j){x2,t), (j?{xi,t)(t?{x2,t) etc. 
We may write 

(0) = jv[<P]0[<P\ jv[(f>o]p[(Po]U*[(P,(Po,t-to]U[(P,(Po,t-to] (1) 

with 

U[<P,<Po,t-to] = y"p[0(r)]e'4^[<^(^)]'^^ (2) 

where the functional integral on the right hand side of Eq. (H) goes between fields (poi^) at 
to and 0(x) at t. The functional p[0o] gives the initial ensemble defining the system. We 
make no assumption of an equilibrium or near-equilibrium ensemble. Later on we shall state 
some general assumptions on p. 
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The separate functional integrals between to and t in U and U* is characteristic of any 



real-time formalism describing the time evolution of a statistical system |]T0|. Rather than 
viewing U*U as determined by two separate functional integrals over the single field (p one 
can introduce two fields, 0_ and 0+, with Lagrangian 



\2 ^ 2 i2 



A 

4!' 



Now 



U*U 



V\(hjV\ 



dr 



(3) 



(4) 



with both 0_ and taking on values </'(x) at t and 0o(^) at to- 

It is convenient to rewrite C in Eq. in terms of new variables and tt defined by |Tl 

1 



^- + 
^2' 



IT 



71 

2 



One easily finds 



(5) 
(6) 

(7) 



where one should be careful not to confuse the in Eqs. with the original field 

appearing in Eqs. (|l])-(|D- In what follows (p "will always denote the variable defined in 
Eq. (^. Noting that tt = at times to and t, one can rewrite Eq. (|l|) as 



{0)= V[<P]0[<P] VIMpIM /P[0(r)]I?[7r(r)]e^^4^[^'-]^- 



with 



L[(p,n] = J £xC[(l),7i:] 



(8) 



(9) 



Our focus is on systems where 0_ and 0+ are large. As we shall see, is naturally large 
while 71 is small, thus one can neglect the 7r^0 term in Eq. (^, which gives 

A 



J^c = df,(f)df,n - m 07r - — 7r0' 



(10) 



But now vr is simply a variable of constraint so that the functional integral over tt in Eq. 
with L replaced by L^, gives 

A 



(□ + m2)0(f,r) + -0=^(f, r) 



(11) 
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which means that only fields satisfying the classical equations of motion 



(12) 



contribute to the functional integral in Eq. (^. In other words, if one defines the evolution 
kernel 



K\ 



bo,t- to 



^IL Lc[<t>,TT]dT 



(13) 



where, as usual, (po and (p are the fields between which the functional integral is taken, then K 
is nonzero only for (p statisfying the field equation with the initial condition (j){tQ, x) = 4>o{x), 
i.e., 

^ ~ i^[0,0o,t-to] = 0, (14) 



{n + m^)<P{x,t) + -<P%x,t) 



and 



lim [0, 00, - ^o] = TT 5[(p{x) - 0o(x)] 



(15) 



The classical nature of K is manifested in Eqs. ([T^) and (15). In terms of K one can write 

(O) as 

(O) = \ P[0] I?[0o] O[0] [0, 00, t - to] p[0o] . (16) 



In fact, the functional integral (|T3D with the Lagrangian of the type (|T0|) provides a conve- 



nient starting point to a diagrammatic approach for classical statistical systems fT^. 



III. EQUATIONS FOR THE GREEN'S FUNCTIONS 

In this section we derive equations for the Green's functions of a system with Lagrangian 
given by Eq. ([To|) . We begin our discussion with the Green's functions (propagators) of the 
free theory corresponding to 

-^^0 = d^cpd^Ti - ni^cp-n , (17) 

then we go on to the equations for the full (classical) theory governed by the Lagrangian (p!0[). 
We limit ourselves in this section to vacuum propagators, leaving the generalization to the 
medium case to Sec. 
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A. The free propagators 



Perhaps the easiest way to get the propagators for the Lagrangian (0) is to go back to 
the 0_ and 0+ variables, in which case 

, 1 



where 



-5 V'; 



and TT are related according to Eqs. (Bf) and (0). Then 



G^^l{x) = (O|T0_(a;)0_(O)|O) 



-ik-x 



(27r)4 F - m2 + ie 
is the usual Feynman propagator. In momentum space 

.(0) . . i 



The propagator for the 0+ field is given by 



p2 _ rfj^l _|_ 



(18) 



(19) 



(20) 



(21) 



with the ie choice reflecting the fact that the 0+ originally come from complex conjugate 
amplitudes. While the Lagrangian ( pISD does not by itself require a or a G\__ propagator, 
the functional integral does naturally give such propagators because of our identification 
0+ = 0_ at the largest time to which we evolve the system. Such propagators also naturally 
occur as Feynman lines which pass from an amplitude to a complex conjugate amplitude in 
the two-time formulation of a calculation of an expectation value or transition probability. 
Clearly, 

Gf_{p) = 27re{po)6{p'-m'), (22a) 
GL°i(p) = 27Te{-po)S{p^ - m^) . (22b) 

Now it is straightforward to go to propagators in terms of and vr fields. For example. 



G 



(0) 
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G 



(0) 



and 



G 



(0) 
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G 



(0) 



(0) 



G 



(0) 



G 



(0) 



etc. One easily finds 

Gf^{p)=nSip' 



m 



G 



G?]ip) = 



(0) 

++ 



G 



(0) 



G 



(0) 



(23a) 



(23b) 



B. Equations for the full propagators 

It will be useful to have a pictorial notation for the various propagators. For the free 
propagators this notation is given in Fig. |l]. For the full propagators our pictorial notation 
is given in Fig. 0. G22 = even in the presence of radiative corrections because of the fact 
that G12 and G21 are causal and anticausal, respectively. (More precisely, one can follow 
a G^2i propagator through a graph much as one can follow an electron through a graph in 
QED. In G22 we can follow Gf^ propagators through the graph and they must form at least 
one closed loop which is impossible because of the anticausal property.) 

p(0)_ ^ 

Cjii - < 

(jj2 — *^VWWVW 

Cj21 — AVWVW^A-* 

FIG. 1: The free propagators 




FIG. 2: The full propagators 



One can write equations for G in terms of one-particle- irreducible parts, S's, which are 
useful in applying perturbation theory and in deriving a Boltzmann equation. For Gu this 
equation takes either the form 

Gn{x,y) = -i dw dz iGf^ {x,w)T,i2{w, z)G2i{z,y) + Gf^ {x,w)T.22iw, z)G2i{z,y) 

10 ^^^^^ 

+Gfj{x,w)J:2i{w,z)Gu{z,y)j + Gf^{x,y) 

or 

Gn{x,y) = -i dw dz \ Gi2{x,w)I:2i{w, z)Gf'-^ {z,y) + Gi2{x,w)Tj22{w, z)G^2i i^^v) 

^ ^^^^^ 

+Gu{x,w)J:i2{w, z)G'S {z,y)j + G'S{x,y) 
7 



+ 



(a) 



+ 1 



+ 



+ © (E 

(b) 



+ 1 



FIG. 3: Diagrammatic represetation of Eq. (25) 



as illustrated in Figs, ^a and Fig. respectively. 

The —i in Eqs. ([25| ) reflects the usual convention that — be the sum of one-particle 
irreducible graphs. Applying + to Eq. ( p5a| ), + to Eq. (|25b ) and substracting, 
one obtains 



(n^; - ay)Gu{x, y) = Jdz {Gu{x, z)J:i2{z, y) - ^21(0;, z)Gii{z, y) 

+Gi2{x, z)J:22{z, y) - S22(x, z)G2i{z, y)} . 

Now it is convenient to write 

d^p f X + y 



(26) 



G{x,y) 



(2vr) 



G 



P e 



-ip-(x-y) 



(27) 



with similar formalae for S(x, y) in terms of E(^-j^,p). We now assume that the dependence 



of G and S on is slow enough that one can replace by x in Eq. (127|). Then using 



x+y 



9 



(28) 



one gets from Eq. (p6D 

■ — G'ii(x,p) = G'ii(a;,p)[S2i(x,p) - Si2(x,p)] + S22(a;,p)[G2i(a;,p) - Gi2{x,p)] . (29) 

It is important to emphasize that Eq. ( p9D is an equation for a classical fleld theory. 
We are dealing with the Lagrangian (|TUp which is classical. In the quantum case, with L 
given by Eq. (0), G22 and En are still zero, but when using the equations derived in this 
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section one must verify that Gu ^ Gu, G21 so that the approximation (j) ^ n which leads 
to the classical theory is valid. This will, of course, depend on the initial ensemble, the 
p in Eq. (|1^), as well as on how far from the initial ensemble the system has evolved. In 
particular we suppose the initial ensemble, p, is dominated by large 0o contributions. Of 
course a system near kinetic equilibrium cannot be described by a classical field theory. 
Turn now to Gu- Analogously to Eqs. one may write 

Gi2{x,y) = —i J dw dz G'f2 {x,w)T.2i{w.z)Gi2iz, x) + G^i2 i^^y) (30a) 

and 

Gi2ix,y) ~ J dw dz Gi2{x, w)'E2i{w .z)G'f2 (z, x) + G^^^ (a;, (30b) 
illustrated in Fig. |. It is now straightforward to get 

2tp-^Gi2{x,p)=0 (31) 

with an identical equation for G2i- Equation (|3T1) does not mean that the graphs shown in 
Fig. ^ do not change Gu from the value given in Eq. (p^ for (^[2^ but it does mean that 
the evolution of the system toward equilibrium is controlled by the equation for Gu where 
the right hand side of Eq. ( ^91) will shortly be identified with the Boltzmann collision term. 

i (^^Gy^vvvs, - ^ (^^G^^VvNA/vv + i - 

(a) 



(b) 

FIG. 4: Diagrammatic represetation of Eq. (pO) 
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IV. EQUIVALENCE BETWEEN THE BOLTZMANN EQUATION AND CLASSI- 
CAL FIELD EQUATIONS 

We now turn to the task of showing that the right hand side of Eq. (|29D is exactly the 
Boltzmann collision term. We do this explicitly only at lowest order in perturbation theory 
for the collision term, after which the higher order corrections should be straightforward. 
We begin by identifying a part of Gu with the Boltzmann phase space particle density. 
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A. The phase space density 



The generalization of Eqs. (pl])- p2BD to a medium is straightforward. The results are, 
in form, exactly like those familiar finite temperature field thoery. One has 

i 



G—ix,p) 
G++{x,p) 



p2 _ ^2 _|_ 

—i 

p2 _ ^2 _ 



G^4x,p) = 27r5{p'-m')[e{-po) + f], 



G-ix,p) = 27T5ip''-m'')[e{po) + f] 



(32a) 
(32b) 
(32c) 
(32d) 



but in contrast to finite temperature field theory / is not the thermal phase space distribution 
but rather 

f = f{x,p,t) (32e) 

is a time-dependent phase space density of particles characterizing the medium. Equa- 
tions (|32D are appropriate to free particles with a phase space density given by /. They are 



not compatible in detail with Eqs. (25) and ([30|). For example, space and momentum de- 
pendent mass effects are not included. Nevertheless, when / ^ 1/A the one-loop correction 
to the mass square is small compared to either the bare mass square or the square of the 
typical particle momentum. Hence we believe Eqs. ( |32D are adequate for our purposes of 
demonstrating the equivalence between classical field theory and the Boltzmann equation 
and match well with the level of accuracy of our Boltzmann equations, Eqs. ( P^D and 
Using equations identical to Eqs. (|23a|) , ( p3b|) etc., one easily finds 

1" 



6-11(2;,^) = 2tt6{p' 
Gi2{x,p) 
G2i{x,p) 



m 



f + 



p2 _ ^2 _|_ ^gp^ 

i 



p2 _ rj^2 _ j^^p^ 



G22{x,p) = 



(33a) 
(33b) 
(33c) 
(33d) 



^ For quantities sensitive to soft modes, like the bulk viscosity, medium corrections to the mass are impor- 
tant 0. 
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and we note that only Gu has changed from Eq. (0). When / is large only Gu is large, 
which means that (j) is large while vr is small, as we have assumed in the discussion before 
Eq. (|I0|). Now that we have found the relationship between the phase space density of 
particles, /, and Gu, we can view Eq. ( pUj ) as an equation for /. Indeed we can rewrite 
Eq. as 

1" 



S21 — S 



12 



/ + 



^22 



(34) 



2 J 2uj{p)' 

where po = ^ip) = a/p^ + m? is to be taken in the S's in Eq. (|3^). Equation (34) has the 



form of a Boltzmann equation. Our task now is to show that the right hand side of Eq. 
agrees with the collision term in 0^ theory, at least when / ^ 1. 



B. The collision term in (/>^ theory 



We are now going to give the lowest order contrubution to the collision term calculated 
directly from the elastic scattering cross section from the graph shown in Fig. |^. The first 
term on the right hand side of Fig. ^ is the gain term for a particle of momentum p while 
the second term is the loss term. With the notation 



^ + V ) /(x,p,t) = G(f,p,t) 



(35) 



one has 



5[uj{p-k)+uj{P)-uj{p)-uj{P-k)] 



{2n)^2uj{p-k)2uj{P)2iu{p)2uj{P-k) 



X 



(36) 



X f{p-k)fiP)[l + f{p)][l + f{P-k)] - f\p)f\P-k)[l + /(p-A;)][l+/(P)] 



C(p) = 





FIG. 5: The lowest order collision term 



The first factor on the right hand side of Eq. (|36D , the |, is a symmetry factor. The term 
quartic in the /'s cancels in Eq. ( pBf ) and one is left with cubic and quadratic terms. The 
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cubic term is 

{ }3 = f{p-k)f{P)[f{p) + f{P-k)] - f{p)f{P-k)[f{p-k) + f{P)] (37) 

with the first two terms on the right hand side of Eq. (^) coming from the gain term and 
the second two terms coming from the loss term. The quadratic term is 

{ }2 = f{p-k)f{P) - f{p)f{P~k) (38) 

We shall see that the classical field evolution reproduces the { part of the collision term 
and, perhaps surprisingly, even the { }2 part. 



C. Evaluating the collision term from classical field theory 



Now we turn to evaluating the right hand side of Eq. (|3^) at lowest order, order A^, 
in perturbation theory. We begin with the S22-term whose lowest order contribution is 
illustrated in Fig. |. One has 

- - {ki + k2 + ks-p}Gu{ki)Gu{k2)Gn{k3) , (39) 



2uj{p) 2uj{p) 3! J (27r)8 
where the second factor of i on the right hand side of Eq. (|39|) comes because ~iT,22 is given 
by the usual Feynman rules, and the 1/3! is the symmetry factor for the graph. In one of 
the three factors of Gu we take /cq < while in the remaining two factors /cq > 0. (This 
is the only way to satisfy the 5-function constraint. The choice of one of the lines to have 



fco < also introduces a counting factor of 3.) Using Eq. ( |33a ), it is now straightforward to 
get 

i^22 ^ /" e-r / i \ i n\ / \ / n i \^ d"^ P d?k 

^ = -y ^[^(p-fc)+^(P)-^(p)-.;(P-A:)] ^^^^^^^^^_^^^^^^^^^^^_^^^^^^^ ^^^^ 

[f{p-k) + \] [f{P) + |] [f{P-k) + |] , 

where we have chosen, say, ki = P, ^2 = P — k, k^ = —{P — k), to match the picture of 

incoming and outgoing lines in Fig. ^. One easily sees that the cubic term in / agrees with 

one of the gain terms in Eq. (p6|). 

Now we turn to the remaining term on the right hand side of Eq. (plj), with S21 — S12 

illustrated in Fig. 0, at order A^. One has 

-i v,i X , 11 /V. V. N •(-^A)2 fd^kid^d^ks ^ ,7,7 a 

Mp) + ^] - = Mp)'^ J i2nr ^ + k2 + ks-p) ^^^^ 

2Txe{k;^)5{kl- m^)G,,{k,)Gn{k2) 
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FIG. 6: Lowest order contribution to 





FIG. 7: S21 — S12 at lowest order 

When > there are two (identical) terms, one having kiQ > and < and the 
other having kiQ < and k2o > 0. When k^Q < it is necessary that kiQ and ^20 both be 
positive. Thus one finds 

' 2 

d^Pd^k r „, s n ^42) 



Mp) 



[f{p) + i] (S21 - S12) 



5[iu{P-k) + iu{P) - uj{p) - iu{P-k)] 

[fip) + I] 



X 



{27t)^2uj{p)2uj{p - k)2uj{P)2uj{P - k) 
X {2 [f{P-k) + i] [/(P) + i] - [f{p-k) + i] [/(P) + \] ] 

where in the first term in { } in Eq. ( ^21) we have taken k^ = p — k, ki = —{P — k), k2 = P 
along with ki ^ ^2, while in the second term we have taken ki = {p—k), k2 = P, and 
^3 = —{P — k) so that the variables match those in Eq. (|3^). We note that the two loss 
terms in Eq. (|37D are in fact identical. 

Now, taking the terms cubic in / in Eqs. (1^) and (0) exactly reproduces the cubic term, 
(^), in Eq. (^). For the quadratic terms we find, from Eqs. (^OD and (0), the result that 
replaces { }2 in Eq- (|38| ) is 

{ }2 = l{fiP-k)f{P) + f{p-k)f{P-k) + f{P)f{P-k) - 2f{p)f{P-k) 

- 2/(p)/(P) - 2f{P-k)f{P) + f{p)f{p-k) + f{p)f{P) + f{p-k)f{P)_ 

Using the fact that f (p—k) f (P—k) and f{P)f{P—k) are equivalent expressions, after inte- 
gration in Eq. (|40|) or (^2]), as are f{p)f{P) and f{p)f{p — k) one finds { }2 = { }2- Finally 



(43) 
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there are the hnear terms in the /'s coming from Eqs. (1^) and (^21) for which there are no 
counterparts in Eq. ( pG]) . In fact the hnear terms do not cancel in Eq. (^) and (^2[) so that 
the classical field theory does not exactly reproduce the Boltzmann collision term. 

Finally, a comment on the level of accuracy at which one can expect the classical field 
theory to reproduce the Boltzmann equation. We always suppose that / 1 but that 
A/ ^ 1 so that a perturbative discussion makes sense. The size of our fields then are 
(p ~ ^// while TT ~ 1 / In going from the full Lagrangian, given in Eq. (0), to the 
classical Lagrangian, given in Eq. (|T0|) , we have dropped a term —^n^cj) compared to the 
interaction term — |y7r0^ which has been kept. Now an estimate of the size of these two 
terms is 

X7r(j)^ ~ A/ (44) 

and 

Xn^cp ~ J (45) 

Thus, the quantum corrections are naturally down by a factor of 1//^ as compared to a 
classical evaluation at a corresponding order of A. Since the maximal possible terms in the 
collision term are of size X^f^, one expects quantum corrections to occur at a level of X^f 
leaving the A^/^ terms still in the classical field theory domain. So in retrospect, our ability 
to get both the X^f^ and X^f^ terms in Eq. ( PBD correctly was to be expected. This means 
that the 2 —>■ 2 parts of a collision term can be obtained from a classical theory calculation. 
This cannot be expected to be the case for 2 — > 4 processes where terms of size A'^/^, A'^/^, 
A^/^, and A'^/^ in the collision terms will occur at leading order in A. One can here expect 
to get the A^/^ and A^/^ terms correctly in a classical field theory, but not the A^/^ and 
A^/^ terms. 
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